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Abstract 

The effect of mixture lengths of vehicles on the asymmetric exclusion model is stud- 
ied using numerical simulations for both open and periodic boundaries in deterministic 
parallel dynamics. The vehicles are filed according to their length, the small cars type 1 
occupy one cell whereas the big ones type 2 takes two. In the case of open boundaries 
two cases are presented. The first case corresponds to a chain with two entries where 
densities are calculated as a function of the injecting rates al and a2 of vehicles type 1 
and type 2 respectively, and the phase diagram (al,a2) is presented for a fixed value of 
the extracting rate (5. In this situation the first order transition from low to high density 
phases occurs at al + a2 = f3 and disappears for a2 > f3. The second case corresponds 
to a chain with one entry, where a is the injecting rate of vehicles independent of their 
nature. Type2 are injected with the conditional probability aa2, where < a2 = na < a 
and n is the concentration of type2. Densities are calculated as a function of the injecting 
rates a, and the phase diagrams (a, (3) are established for different values of n. In this 
situation the gap which is a characteristic of the first order transition vanishes with in- 
creasing a for n / 0. However, the first order transition between high and low densities 
exhibit an end point above which the global density undergoes a continuous passage. The 
end point coordinate depends strongly on the value of n. In the periodic boundaries case, 
the presence of vehicles type2 in the chain leads to a modification in the fundamental 
diagram (current, density). Indeed, the maximal current value decreases with increasing 
the concentration of vehicles type 2, and occurs at higher values of the global density. 
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1 Introduction 



Without doubt, an efficient transportation system is essential for the functioning and 
success of modern industrialized societies. But the days when freeways were free are 
over. The increasing problems of roadway traffic raise the following questions: Is it still 
affordable and publicly acceptable to expand the infrastructure? Will drivers still buy 
cars when streets are effectively turned into parking lots? Automobile companies worried 
about their future market, have spent considerable amounts of money for research on 
traffic flow and on how the available infrastructure could be used more efficiently by new 
technologies. Indeed, traffic flow is an interesting field of inter-disciplinary research; it has 
attracted the interest of many researchers from different disciplines. Like mathematicians, 
chemists, and engineers, physicists have addressed the problems of traffic flow many years 
ago and they have been trying to understand the fundamental principals governing the 
flow of vehicular traffic using theoretical approaches based on concepts and techniques of 
statistical physicsfl, 2, 3]. The approach of a physicist is usually quite different from that 
of a traffic engineer. A physicist tries to develop a model of the traffic by incorporating 
only the most essential ingredients which are absolutely necessary to describe the general 
features of typical real traffic. There are two different ways for modelling traffic: The 
macroscopic models which are based on fluid-dynamical description and the microscopic 
ones where attention is explicitly focused on individual vehicles which are represented 
by particles. The interaction is determined by the way the vehicles influence each others 
movement. In other words in the microscopic theories traffic flow is considered as a system 
of interacting particles driven far from equilibrium. Thus, it offers the possibility to study 
various fundamental aspects of the dynamics of truly non equilibrium systems which are 
of current interest in statistical physics [4, 5]. Within the conceptual framework of the 
microscopic approach, the particle hopping models describe traffic in terms of stochastic 
dynamics of individual vehicles which are usually formulated using the language of cellular 
automata (CA)[6]. In general, CA are an idealization of physical systems in which both 
space and time are assumed to be discrete and each of the interacting units can have 
only a finite number of discrete states, thus in CA models of traffic the position, speed, 
acceleration as well as time are treated as discrete variables and the lane is represented by 
a one-dimensional lattice, each site represents a cell which can be either empty or occupied 
by at most one vehicle at a given instant of time. The computational efficiency of CA is 
the main advantage of this approach. Much effort has been concentrated on stochastic 
CA models of traffic flow first proposed by Nagel and Schreckenberg[7]and subsequently 
studied by other authors using a variety of techniques [8, 9]. The stochastic dynamics of 
interacting particles have been studied in the mathematical and physical literature [10]. In 
the physical case, driven lattice gases with hard core repulsion provide models for diffusion 



2 



of particles through narrow pores and for hopping conductivity [11], and belong to the 
general class of non-equilibrium models which includes driven diffusing systems [12, 13]. 
They are closely linked to growth process[14, 15], and can also be formulated as traffic jam 
or queuing problems [16]. The fully asymmetric exclusion model (FAEM) which is also 
described in literature as the Totally asymmetric exclusion process (TASEP) corresponds 
to the case where particles hop only in one direction. This model can be divided into 
four classes according to the dynamics (sequential or parallel) and the choice of boundary 
conditions (open or periodic). In the sequential dynamics in which each particle has a 
probability At(time interval) of jumping to its right-hand neighbour if this neighbouring 
site is empty, has been solved exactly in one dimension with open boundary conditions[17, 
18]. Recently Shaw et al [19] studied the protein synthesis using the TASEP with extended 
objects i.e. a system with particles of length I > 1 using numerical and analytical tools for 
random sequential updating. The stochastic parallel update of the asymmetric exclusion 
model was studied using Monte-Carlo simulations [20], however its deterministic case i.e. 
particles moving forward with probability q = 1, was studied analytically by Tilstra and 
Ernst[21]. The aim of this paper is to study the effect of mixture lengths I = 1 and / = 2 
of vehicles on the FAEM in the case of deterministic parallel dynamics in both open and 
periodic boundary conditions using numerical simulations. This paper is organized as 
follows; section 2 is devoted to explain the model with three cases i.e. open boundaries 
with two entries, open boundaries with one entry and periodic boundaries. In the third 
section we present the main results obtained in these cases with a critical discussion. 
Section 4 is devoted to a conclusion. 

2 Model 

We consider a chain of N sites in which the vehicles are filed according to their length. 
However, vehicles of typel means short cars (I = 1) which occupy one site and type2 the 
long vehicles (I = 2) occupy two sites. Type2 moves by one site at each time step even if 
two sites ahead are empty, the same applies to typel. Hence Fig.l displays an example 
of some configurations. 

2.1 Open boundaries 
2.1.1 The case of two entries 

In this case, typel are injected in the first site (the first entry) with an injecting rate al, 
while type2 are injected in the third site (second entry) with an injecting rate a2, if these 
entries are empty. However, in order to avoid the on-ramp at the moment when type2 are 
injected, we should take into account the following conditions: 
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The second site must be, also, empty, to not hinder the movement of typel present in the 
first site at the next time step, because type2 enters the road on two time steps due to its 
length. Thus, the distance between the two entries should be at least equal to 1. 

2.1.2 The case of one entry 

For both practical and theoretical reasons, some times different boundary conditions are 
used. Imagine a situation where a multilane road, containing different kinds of vehicles 
(car and trucks), is reduced to one lane due, for example, to a road construction. Such 
a situation can be modelled by using our model. In this case the multilane part of the 
road acts as a particle reservoir (both types), and the one part lane is represented by a 
chain with N sites. Here the vehicles will be injected at the same entry depending on 
the injecting rate a of vehicles and the concentration n of vehicles type2. In this case, 
numerical simulations are limited only when < n < 0.5, since in the road we have a few 
type2 vehicles in comparison to typel (we assume that type2 are trucks and typel are 
cars). Hence, we denote by q2= n«, the injecting rate of vehicles type2. However, the 
process of the injection of the vehicles is as follows: 

If the first site is empty, a random number < R < 1 is chosen. Then, if R < a, a vehicle 
is injected. But its nature depends on the value of a2; hence, if q2 < R < a, then typel 
could enter the road; while, if R < a2, type2 can enter the road. 
Note that < q2 < a, in contrast with the first boundary case where < a2 < 1. 

2.2 periodic boundaries 

The periodic boundary conditions for the FAEM with one species of particle on a ring has 
been studied exactly by Schadschneider and Schreckenberg[8, 9]. Here we study the FAEM 
with two types of vehicles for periodic boundaries with deterministic parallel dynamics. 
To illustrate this situation let us consider a ring of N sites where C\ and C2 are the 
densities of sites occupied by the vehicles typel and type2 respectively. Hence, the global 
density is given by C = C\ + C2. We assume that n = is the concentration of the sites 
occupied by type2. However, C\ = (1 — n)C and C2 = nC. The main idea of studying 
this case is to establish the fundamental diagram (I,C) for several values of n. The mean 
value of the current in the steady state is given by: 



In our computational studies we have considered a chain with N = 1000 sites. Simulation 
of the closed system begins with particles randomly distributed around the ring depending 
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3 Simulations and results 
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on their densities C\ and C2, whereas the open system begins with a small number of 
vehicles randomly distributed in the chain. The systems run for 20 000 MCS to ensure 
that steady state is reached for the periodic case and 40 000 MCS for open systems, at 
this moment data including the current and density are collected. In order to eliminate 
the fluctuations 25 initial configurations were randomly chosen. 

3.1 Open boundaries 

3.1.1 The case of two entries 

We recall that our aim in this case is to explore the phase diagram (al,a2) for a fixed 
values of the extracting rate (3. Fig.2 shows the profile of the global density p (which 
is the density of all sites occupied by both types of vehicles) versus al for (3 = 0.3 and 
various injecting rates a2. We distinguish two regions: 

The first region is when (al + a2) < (3 in which p increases with al, but this increase 
depends on the value of a2. Indeed: 

p{al, a2 = 0) < p(al, a2 = 0.1) < p(al, o2 = 0.2) 

On the other hand this region corresponds to the low density phase due to the low values 
of p. The second region occurs when (al + a2) > (3 in which p become constant (about 

0. 77) and does not depend on al and a2. This is the high density phase. 

At (al + a2) c = (3 the global density p is discontinuous and undergoes a jam between the 
two regions which is a characteristic of the first order transition. 

We can see that the system undergoes the usual first order transition at al = (3 for 
a2 = 0.0 in agreement with[22]. Also we remark that the density takes its higher values 
when a2 > (3 even for al = 0.0, this means that there is no low density phase in this case 

1. e. the system prefers to be in the high density phase. In fact typel and type2 are not 
correlated at the beginning of the road due to the conditions mentioned in section (2.1.1), 
this explains why the first order transition occurs at(al + a2) c = (3. For a2 > (3 this first 
order transition disappears due to the greater proportion of type2 in the road, indeed the 
number of type2 which enter the chain is greater than those which leave it leading to high 
density phase independently of al. Collecting all of these results we obtain the phase 
diagram (al, a2) for (3 = 0.3 in Fig.3 The same results can be obtained for others values 
of (3. 

3.1.2 The case of one entry 

This situation corresponds to the case when typel and type2 enter the chain in the same 
entry. We recall that the FAEM with open boundary conditions and deterministic parallel 
dynamics exhibits a first order transition [22], in which p is discontinuous and undergoes 
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a gap between the low and high density phases at a c = (3. In our model the gap decreases 
with increasing a and vanishes for its higher values. Fig. 4 shows the profile of the global 
density p versus a for a2 = na and various extracting rates /3, here n=l/4 corresponding 
to the concentration of type2. We remark that for (3 = 0.0 — 0.3 the first order transition 
occurs at a c = (3 in agreement with [22] and at a c < [3 when (3 = 0.4 — 0.7. For 
(3 = 0.8 — 0.9, there is no phase transition. In fact, in this case about 1/4 of the vehicles 
which enter the road could be type2, and their number becomes important as long as 
a increases, we recall that types2 occupy two sites due to their length so, it is clear in 
this situation that the global density in the road becomes larger than the n=0 case (i.e. 
the chain contains only typel vehicles). Such a comparison is illustrated in Fig.5, in 
which the n=l/4 density gap is smaller than the n=0 one. The first order transition 
occurs ata c < /3(n=l/4) instead of a c = /3(n=0.0). This result is not observed when 
(3 = 0.1 — 0.3(Fig.4), because the system reach the high density phase for small values of 
q, which means that the probability for type2 to enter the road is small. Furthermore, 
for (3 = 0.8 — 0.9, the gap vanishes (Fig.4) and we have a continuous passage from the low 
to the high density phase. In Fig.6 we display the (a, (3) phase diagram for three values 
of n. It is found that the first order line transitions, separating the low and high density 
phases in the case of n = 0.0 [22], is modified when n = 0.25, 0.5. Indeed, for n = 0.25 and 
n = 0.5 the line transitions terminate by end points above which the first order transitions 
disappear. These end points are located at (a = 0.7,(3 = 0.74) and (a = 0.5, (3 = 0.57) 
for n = 0.25 and n = 0.5 respectively. 

3.2 Periodic boundaries 

The periodic boundary conditions for the FAEM with one species of particle (Typel) on 
a ring has been studied exactly by Schadschneider and Schreckenberg[8, 9], where the 
relation between current and density is given as follows: 

I(C,P) = 1/2[1 - y/l-4qC(l-C)] 

p is the braking probability and q = 1 — p is the hopping rate. 

An interesting feature of this expression is that the current is invariant under the oper- 
ation C — ► (1 — C) which interchanges particles and holes. Therefore the fundamental 
diagram is symmetric about C m = 1/2. This symmetry is conserved for all p, but breaks 
down for a fixed value of p and V max > 1 where the magnitude of C m decreases and the 
current increases with increasing V max . 

In our model, introducing a small concentration of type2 leads to a breaking of the sym- 
metry between the holes and particles and the expression of the current I(C, P) mentioned 
above doesn't work. Fig. 7 shows the fundamental diagram (I,C) for different concentra- 
tions n of type2, it is clear that the magnitude of C m increases and the current decreases 
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when n increases. We can say that so long as C is sufficiently small the vehicles are too 
far apart to interact mutually, therefore the current increases with C but this increase 
is controlled by the concentration n. Indeed type2 takes lot of time in passing the local 
detector due to its length in contrast to typel, which reduces the current and shifts to 
the right C m . In the particular case n = 1, with a random sequential dynamic, Shaw et 
al[19] found using analytical and numerical calculations, that the current I isn't symmet- 
ric about C m = 1/2. Indeed, the maximal current is lowered from 0.25(n=0, the chain 
contains only typel vehicles) to I m = 0.1715(n=l) and the density C m is shifted from 
0.5 to C m = 0.585. However, using numerical simulations with a deterministic parallel 
dynamics update, we find that the maximal current is lowered from 0.5(n=0)[8, 9] to 
I m = 0.332 (n=l) and the density C m is shifted from 0.5 to C m = 0.669 as presented in 
Fig.7. 

We note that Fig. 7 can be reproduced analytically from the expression of the current 
cited above (for q=l) by rescaling the density C and the current I. Indeed, instead of C 
one then has C(l — n/2)/(l — nC/2) and for the current I one gets 1/(1 — nC/2), hence 
the expression of the current becomes 1(C) = 1/2[1 — nC/2— \ 2C — 1 — nC/2 |]. 

4 Conclusion 

We have studied the effect of a mixture of lengths of vehicles on the traffic flow using 
numerical simulation with a deterministic parallel update. For this purpose we have 
defined two types of vehicles depending on their length, typel the small cars (Z = 1) 
which occupy one site, and type2 the long ones (I = 2) taking two sites. We have 
investigated both open and periodic boundaries. In the former case, two different cases are 
considered. In the case of two entries, the global density p is discontinuous and undergoes 
a jam between the low and high density phases which is a characteristic of the first order 
transition at (al + a2) c = (3; and for a2 > (3 the system is in the high density phase 
then no phase transitions is observed. In the case of one entry, the concentration of type2 
plays a major role, indeed the gap between low and high density phases decreases when 
increasing a and vanishes for higher values. Furthermore the phase diagram exhibits an 
end point above which the first order transitions disappear. 

In the periodic boundary case, the presence of type2 vehicles leads to a breaking of the 
symmetry holes-particles, then the magnitude of c m increases, shifting to high density 
and the maximal current decreases. 
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Figure captions: 

Fig. 1: Example of configurations obtained after three steps for system size L = 8 
where both typel and type2 move with V max = 1. 

Fig. 2: The variation of the global density p versus al for different values of a2 and 
(3 = 0.3 in the case of open boundaries with two entries, the number accompanying each 
curve denotes the values of a2. 

Fig. 3: The phase diagram (al,a2) for (3 = 0.3 in the case of open boundaries with 
two entries. 

Fig. 4: The variation of the global density p as a function of a with a2 = l/4a 
and for different values of (3 in the case of open boundaries with one entry. The number 
accompanying each curve denotes the value of (3. 

Fig. 5: The variation of the global density p as a function of a with a2 = 1 /4a (circle 
dot) and a2 = 0.0a(square dot) for (3 = 0.7 in the case of open boundaries with one entry. 

Fig. 6: Phase diagram in the (a, (3) plane for several values of n (concentration of 
type2). Square and up triangle dots correspond respectively to a2 = 0.5a and a2 = 
0.25a where the lines of the first order transitions ended at the end points (0.5,0.57) and 
(0.7,0.74). Circle dots recall the case a2 = 0.0a. 

Fig. 7: The fundamental diagram (I,C) for several values of n. The symbols are 
obtained by simulation whereas the lines represent the analytical expression of the current. 
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